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Abstract 

Motivated by future collider proposals that aim to measure the Higgs properties 
precisely, we study the partial decay widths of the lightest Higgs boson in the minimal 
supersymmetric standard model with an emphasis on the parameter region where all 
superparticles and heavy Higgs bosons are not accessible at the LHC. Taking account 
of phenomenological constraints such as the Higgs mass, flavor constraints, vacuum 
stability, and perturbativity of coupling constants up to the grand unification scale, 
we discuss how large the deviations of the partial decay widths from the standard 
model predictions can be. These constraints exclude large fraction of the parameter 
region where the Higgs widths show significant deviation from the standard model 
predictions. Nevertheless, even if superparticles and the heavy Higgses are out of the 
reach of 14TeV LHC, the deviation may be large enough to be observed at future e^e“ 
collider experiments. 


1 Introduction 


The discovery of the Higgs boson at ATLAS and CMS experiments DEa made a revolntionary 
impact on the held of particle physics. It not only conhrmed the so-called Higgs mechanism 
for the electroweak symmetry breaking, but also opened a new possibility to perform a 
precise test of the standard model (SM) by studying the properties of the Higgs boson. In 
the SM, the coupling constants of the Higgs boson with other particles are well understood 
using the fact that the masses of quarks, leptons, and weak bosons originate in the vacuum 
expectation value (VEV) of the Higgs held, resulting in the prediction of the partial decay 
widths of the Higgs boson into various particles. 

In models with physics beyond the SM (BSM), measurements of the Higgs couplings 
provide even exciting possibilities. In large class of BSM models, there exist new particles at 
the electroweak to TeV scale, which ahect the properties of the Higgs boson. Thus, with the 
detailed study of the Higgs properties at collider experiments, we have a chance to observe 
a signal of BSM physics. Such a study will be one of the major subjects in forthcoming 
collider experiments, i.e., the LHC and future colliders like ILC and TLEP [3]. 

Low energy supersymmetry (SUSY) is a well-motivated candidate of BSM physics. Com¬ 
pared to the SM, the particle content is enlarged in SUSY models. Even in the minimal setup, 
i.e., in the minimal SUSY SM (MSSM), there exist two Higgs doublets, Hu and as well 
as superparticles. The lightest Higgs boson h, which plays the role of the “Higgs boson” 
discovered by ATLAS and CMS, is a linear combination of the neutral components of Hu and 
Hd, while there exist other heavier Higgses. In the case where the mass scales of the heavier 
Higgses and the superparticles are high enough, the properties of h are close to those of 
the SM Higgs boson. On the contrary, if the heavier Higgses or superparticles are relatively 
light, deviations of the Higgs properties from the SM predictions may be observed by future 
collider experiments. With the precise measurement of the partial decay widths (or branch¬ 
ing ratios) of the Higgs boson, information about the heavy Higgses and/or superparticles 
may be obtained even if those heavy particles can not be directly discovered. 

In this paper, we discuss how low the mass scales of the heavier Higgs bosons and su¬ 
perparticles should be to observe a deviation. We evaluate the partial decay widths of 
the lightest Higgs boson in the MSSM, taking account of the following phenomenological 
constraints: Higgs mass, flavor constraints of the B mesons, stability of the electroweak 
(SM-like) vacuum against the transition to charge and color breaking (CCB) vacua, and 
perturbativity of coupling constants up to a high scale. These constraints exclude large 
fraction of the parameter region giving rise to a significant deviation. Even so, we will see 
that the deviations of the partial widths from the SM predictions can be of 0(1) % for some 
of the decay modes in the parameter region allowed by the above-mentioned constraints. In 
particular, the deviations may be large enough to be observed by future future e+e“ colliders 
like ILC and TLEP even if superparticles are so heavy that they would not be observed at 
the LHC. 

The organization of this paper is as follows. In Sec. O we briefly overview the properties 
of the Higgs bosons in the MSSM. We also summarize the phenomenological constraints that 
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are taken into account in our analysis. Then, in Sec. El we calculate the partial decay widths 
of the lightest Higgs boson in the MSSM and discuss how large the deviation from the SM 
prediction can be. Sec. IHis devoted for conclusions and discussion. 


2 MSSM: Brief Overview 


2.1 Higgs sector of the MSSM 

We review some of the important properties of the Higgs sector in the MSSM. There are 
two Higgs doublets, and H^- As the neutral components acquire VEVs, the electroweak 
symmetry breaking (EWSB) occurs. The ratio of the two Higgs VEVs is parameterized by 
tan/3 = Assuming no CP violation in the Higgs potential, the mass eigenstates 

are classified as lighter and heavier CP-even Higgs bosons (denoted as h and H, respectively), 
CP-odd (pseudo-scalar) Higgs A, and charged Higgs In the following, we concentrate on 
the case where the masses of the heavier Higgses (il. A, and H^) are much larger than the 
electroweak scale. Then, the lightest Higgs boson h should be identified as the one observed 
by the LHC. On the other hand, the masses of the heavier Higgses are almost degenerate. 
We parameterize the heavier Higgs masses using the pseudo-scalar mass m^. 

At the tree level, the lightest Higgs mass is predicted to be smaller than the Z-boson 
mass, while it is significantly pushed up by radiative corrections la El El [71E]. The mass 
matrix of the neutral CP-even Higgs bosons is denoted as 

. .2 _ m| cos^/3- 1 -sin^/3-|-— (m|-|-m^) cos/3 sin/3-|- 
^ — (m|-I-m^) cos/3sin/3-I-m| sin^/3-|-cos^/3-f 5M|2 ’ ’ 


where 5Mfj represents radiative corrections. 

At the one-loop level, the top-stop contribution dominates the radiative correction to the 
lightest Higgs mass, and is approximated as 


Sml ~ 


3mj 


log 


rxir 


0 . ' 9 

mf m~ 

^ t 



( 2 . 2 ) 


where v ~ 246 GeV is the SM Higgs VEV, rrit is the top-quark mass, m? = (with 

and mi^ being the lighter and heavier stop masses, respectively), and Xt = A^ —/r cot /3 (with 
At and /i being the tri-linear scalar couplings for stop and the SUSY invariant Higgsino mass 
parameter, respectively)]^ The top-stop contribution can significantly enhance the lightest 
Higgs mass. On the other hand, the bottom-sbottom contribution to the lightest Higgs 
mass becomes sizable when the bottom Yukawa coupling is large. It is likely to decrease the 
lightest Higgs mass. 

When stop masses are 0{1) TeV, there are up to four solutions for At to satisfy the 
observed value of the Higgs mass rrih, for which we use rrih = 125.7 GeV HD]. Let us call 
these four solutions as 

In this paper, we adopt the convention of the SLHA format [9]. 
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• NS: negative At with smaller \At\, 

• NL: negative At with larger \At\, 

• PS: positive At with smaller \At\, 

• PL: positive At with larger \At\. 

Assuming universal sfermion masses at the SUSY scale, the value of \At\ is typically a few 
times larger than the stop mass for NL and PL cases. Such a large value of |At| has signihcant 
phenomenological implications, as we will discuss in the next section. 

Since the one-loop correction to the Higgs mass is comparable to the tree-level value, 
higher order corrections are necessary to obtain reliable results. In particular, QCD correc¬ 
tion, which appears at the two-loop level, and a large hierarchy between the SUSY scale and 
the electroweak scale require the resummation of the leading and sub-leading logarithms. 
We use FeynHiggs 2.10. 2 dH E2l DSl [a [15] for the precise evaluation of the Higgs masses 
(as well as the mixing parameters and the partial decay widths of h). 

At the tree level, Hu {Hd) couples only to up-type quarks (down-type quarks as well as 
leptons). However, this is not the case once radiative corrections due to superparticles are 
taken into account. The Higgs couplings to bottom quark and tan lepton can be subject to 
sizable corrections even when SUSY breaking scale is very large. Let us parameterize the 
effective hbb and hit vertices including radiative corrections as [m El isi nsiEni ED 

-C,s = Vb e^jbnWdQl + + yt et^inQ^Hi + ^yt inQ^H^/ + h.c., (2.3) 


where b^, tji, and Ql are right-handed bottom, right-handed top, and third-generation 
quark-doublets, respectively. In addition, i, j and k are SU{2)l indices, while the color in¬ 
dices are omitted for simplicity. Here, /Ayt, and Ayt are non-holomorphic radiative corrections 
to the Yukawa coupling constants 


#2 


After the electroweak symmetry breaking, the Yukawa couplings are related to the quark 


masses as 


mb = -^v cos /5 ( 1 H- — tan /? ) = -^v cos 13(1 + Ab), 

V2 \ Vb J \/2 

mt = -^nsin /9 (1 H- — cot /5 ) = -^v sin/5(l -|- A*), 

\/2 \ yt J V2 


(2.4) 

(2.5) 


where, at the leading order in the mass-insertion approximation, Aj is given by 


Ab ~ 


/(m? ,m|,/i^) | tan/?. 


Svr 


( 2 . 6 ) 


For a detailed treatment of the non-holomorphic corrections, see Refs. [221123) . 

These relations hold at the SUSY breaking scale. Thus, the quark masses and Yukawa couplings in the 
formula should be understood as the running parameters at the scale. 
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A* ~ 


Ms/i I {ml , ml , M|) + I (rnl , ml , /x^) 


Svr 


cot /9, 


(2.7) 


with M 3 being the gluino mass. The loop integral is defined as 


I {a, b, c) 


ab In a/b + 6 cln 6 /c + ca In c/a 
(a — b)(b — c)(a — c) 


(2.8) 


Notice that A;, is enhanced when tan (3 is large, while Aj is suppressed by cot [3. 

The mass eigenstates of the Higgs bosons are given by linear combinations of Hu and H^. 
The CP-even parts of their neutral components are related to the mass eigenstates as 


Re(»2) 

Re(//») 



{v sm(3 + h cos a -|- M sin a), 
{v cos {3 — h sin a + H cos a). 


(2.9) 

( 2 . 10 ) 


The mixing angle a depends on the pseudo-scalar mass tua, and shows a decoupling be¬ 
haviour, i.e., cos{(3 — a) —)■ 0 as m^ —)■ cxo. In this limit, h behaves as the SM Higgs boson. 
Using Eq. fl2.ip . we obtain [21] 


cos{(3 — a) 


sin4/j / 6Ml - _ 6Ml \ ( 

2m\ \ 2m|cos2/? m|sin2 /?/ V^a/' 


( 2 . 11 ) 


In Fig. m we show the behavior of cos(/9 — a) as a function of rrxA with the masses of 
superparticles being fixed. Here, all the sfermion mass parameters, rug, m^, mj^, mi, 
and rUg, are taken to be universal at the SUSY scale Mgusv, where these mass parame¬ 
ters are soft SUSY breaking masses of sfermions with gauge quantum numbers of (3, 2,1), 
(3,1,—|), (3,1, i), (1,2,—1), and (1,1,1), respectively, with the numbers in the paren¬ 
thesis being quantum numbers for SU{3)c, SU{2 )l, and f/(l)y. Throughout our study, we 
take MgusY = ■ In addition, the sfermion masses are assumed to be universal 

in generation indices. We can see that radiative corrections can enhance cos(/5 — a) by an 
order of magnitude when is large, while it is comparable to the tree-level value with the 
smaller \At\ solutions. 

Denoting hff coupling (with / being the SM fermions) as 


'^hff = 9hffhff, 


we obtain the hbb coupling constant as 

1 - 


9 hbb 


sm a 


cos (3 
sin(/9 — a) 


- Ab cot a cot (3 (SM) 
1 -F Ab 

tan (3 — Ab cot (3 


1 + Ab 


cos(/5 


a 


„(SM) 

9hbb ’ 


( 2 . 12 ) 


(2.13) 
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Figure 1: cos(/3 — a) is plotted as a function of with tan/? = 35, M 3 = —/i = 5TeV, 
and the approximate GUT relation among gaugino masses (see Eq. fl3.2p i. All the sfermion 
masses are assumed to be universal at Msusy, and are also taken to be 5TeV. The red and 
blue lines correspond to the cases of the PL and PS solutions of At, respectively. The black 
dashed line is the tree-level value. 

where the superscript “(SM)” is used for the SM prediction. When rriA is relatively large, 
sin(/? —a) is almost unity, while the second term proportional to cos(/? —a) induces a sizable 
deviation from the SM value. Similar relation holds for hfr vertex, and A,- is approximately 
given by 


- -^^ 2 htan/?/(m?^,M 2 ^/i^), (2.14) 

with M 2 being the Wino mass. Quantitatively, A,- is smaller than A^ in the parameter space 
of our study. As we will see later, P(h —?■ bb) and P(h —)■ fr) may show sizable deviations 
from the SM predictions even if is above TeV. 

The hit coupling constant is obtained as 


9htt 


cos a 
sin /? 


1 — At tan a tan /? (sm) 
1 +At 


■ cot/? - Attan/? ' 

smi/? — a) H-:-cosi/? — a) 

vv t 1 + At vv t 


(SM) 

9 hit • 


(2.15) 


The deviation from the SM value mainly comes from the second term in the bracket and is 
not enhanced by tan/?, since At is proportional to cot/?. 

The gauge-boson hnal states are also important. In order to calculate the partial decay 
widths of the Higgs to gauge bosons, we have modihed FeynHiggs 2.10.2 package to properly 
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#4 


We 


take account of the effect of non-holomorphic correction to the Higgs interaction, 
have also modihed the package to take the Oefr approximation [2S] for calculating the partial 


widths, in which the renormalization scale of the Higgs wave functions is set to be p'^ = 0, and 
the effects of radiative corrections are included in the mixing between the light and heavy 
Higgs bosons. Then, as we will see below, the partial decay widths show proper decoupling 
behavior in the large rriA limit, 


#5 


The processes induced by triangle loops, h ^ gg and 77 , have been important for the 
Higgs discovery at the LHC. They are also important in studying new particles that couple 
to the Higgs boson. In SUSY models, the stop and sbottom loops contribute to the hgg 
coupling. It is expressed by an approximate formula (cf. Refs. [25l [26l |27]) 


9hgg 

^(SM) 

9hgg 


Qhtt 

(SM) 

9htt 




m. 




+ 




2 ' 2 

m- m. 


2 2 

m - m- 

/l J2 . 


(2.16) 


/■ /= 

up to D-term and bottom-loop contributions. Here, — /itan/3. In the right-hand 

side, the hrst term comes from the top loop, while the second term is given by the stop and 
sbottom loops. The correction is positive in the non-mixing limit (i.e., Xf —)■ 0), whereas it 
becomes negative when the mixing terms are sizable. We also note here that hVV couplings 
(with VV = W~^W~ and ZZ) are approximately given by ghvv/d^^ = sin(/5 — a). This 
ratio is very close to unity, and hence the deviations in these modes are very small. 


2.2 Constraints 

Before discussing the possibility of observing a deviation of the Higgs partial widths from 
the SM prediction at future colliders, we summarize phenomenological constraints on the 
MSSM parameter space which are adopted in our analysis. 


2 . 2.1 Bg —^ pAji 


In the SM, the flavor-changing decay, Bg —)■ , proceeds by virtual exchanges of the Z and 

W bosons. They are suppressed by the hnal-state helicity. In contrast, SUSY contributions 
can be enhanced considerably by large tan/?, when virtual exchanges of the heavy Higgs 
boson contribute to the decay [aniEi]. The branching ratio is expressed as [321 Eal El ES] 


Br(i?s —t 


GW t2 




X 


fBWBjB. WA 


ts \ 




m 


Bs 


Ami 

1 - A 


m 


Bs 


Cqi - a 


Qi| 


+ 


2m, 


(c«2 - 7^) + ^ (c„ - c;„) 

Bs 


(2.17) 


We use Eq. (12.131) for the hbb coupling to calculate the partial decay widths of —>■ 77 , 'yZ and gg. We 
found that the factor of (1 — Af, cot a cot/3) was missing in FeynHiggs 2.10.2 package (see Eq. (I2.13l) i. 

In particular, the partial width T{h ^ gg) converges to the SM value in the limit of large itia and 
squark masses. This behavior looks inconsistent with the result shown in Ref. [29]. 
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where Gp is the Fermi constant, is the i^^-meson mass, fp, is the decay constant of 
Bg, is the muon mass, is the lifetime of Bg, and Vij is Cabbino-Kobayashi-Maskawa 
matrix element. In the above expression, Cio, Cqi and Cq 2 are the Wilson coefficients of the 
effective operators, Oi oc (s 7 ^Pl 6 )(£ 7 ^ 75 £), {sPjib){£i), and {sPRb){i'j 5 i), respectively, while 
C' are obtained by flipping chiralities, R L. Among them, the SM contribution appears 
only in Ciq. Including higher order contributions and taking account of effects of the Bg-Bg 
oscillation, the SM prediction becomes [M] 

Br(P, ^ /iV)sM = (3.65 ± 0.23) x 10"^. (2.18) 

This can be compared with the LHCb measurements 1371 

Br(Ps —)■ )LHCb = (2.8lo;6) ^ (2.19) 

Defining ABr(Ps = Br(i?s —)■ — Br(Ps —)■ /i^/i“)sM, where the first term in 

the right-hand side includes both the SUSY and SM contributions, the 95% C.L. bound is 
estimated as 


-2.3 X 10"^ < ABr(P, ^ /rV) < 0.6 x 10“®. 


( 2 . 20 ) 


We will adopt this constraint in our numerical analysis. 

Even when superparticles are heavy, they affect the branching ratio through non-holomorphic 
contributions to the heavy Higgs couplings. Including radiative corrections and diagonaliz¬ 
ing the quark mass matrices, effective couplings of the heavy Higgs bosons to the down-type 
fermions become [301 l3T] 




grrib 


A 


FC 


\/2mw cos /3 (1 + l^ft) (1 + l^o) 
grrii 1 


VtbVZisLbR) (H + iA) 


\/2mw cos /3 1 + A^ 


{^l^r) {H + iA) + h.c.. 


( 2 . 21 ) 


where Aq = A;, — Apc- The flavor-changing coupling is induced by Apc as 


Afc = Yl^/iAitan/? J(m|,m|,/i^). (2.22) 

Here, soft scalar masses are assumed to be universal in generation, and A^ is obtained by 
substituting f —)■ £ in Eq. fl2.14p . Then, the Wilson coefficients receive Higgs-mediated 
contributions. 


c'qi ^ -a 


Q2 — 




tan^ {3 fiAi 


9wmw^\ (1 + 


^iii^^ini 1 )y 


(2.23) 


where rrib is the bottom-quark mass, 9w is the Weinberg angle, and xi^ = The 

non-holomorphic correction Apc as well as Aj, does not decouple even for very heavy super¬ 
particles. We will see later that the corrections to Br(Ps —^ 9'~) can be sizable and that 
the constraint excludes some part of the parameter space of our interest. 
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The branching ratio of the inclnsive decay of 6 —)■ sy may also be sensitive to the non- 
deconpling contribntions to the (charged) Higgs boson. In the numerical analysis, dehning 
ABr(6 —)■ sy) = Br(6 —)■ sy) — Br(6 —)■ sy)sM, we adopt the 95% C.L. bonnd, 

-3.6 X 10"® < ABr(6 ^ sy) < 9.2 x 10“^ (2.24) 

where the experimental valne Br(6 —)■ sy)exp = (3.43 ± 0.21 ± 0.07) x lO”"^ [S^ and the SM 
prediction Br(6 —)■ sy)sM = (3.15±0.23) x 10“^ [32] are combined. At the cnrrent accnracies, 
Bg —^ imposes more stringent bonnd on the parameter space than b ^ sj except when 

snperparticles are light. 

In the nnmerical analysis, Superlso 3.4 HDl is nsed for evalnating the SUSY contribntions 
to the branching ratios as well as the SM predictions. In onr analysis, we assnme that the 
sqnark masses are nniversal in generations. If the sqnark masses are non-nniversal, there 
are extra contribntions to Apc, and the flavor constraints are affected (see e.g.. Ref. m)- 
Snch a non-nniversality is expected even in the model with the nniversal scalar masses at 
the GUT scale. Thns, it shonld be noted that the flavor constraints that we will show below 
are jnst for a particnlar choice of the sqnark-mass parameters and may change if the sqnark 
mass matrices have non-nniversal structnres. 


2.2.2 Vacuum stability 


With sufficiently large \At\, CCB vacua arise, and the minimum of the scalar potential with 
the correct EWSB becomes a false vacuum [321 SSI SH SSj- When |p| <C |Ai|, stop and the 
up-type Higgs helds acquire large VEVs at the CCB vacua, while the VEVs of other helds 
are relatively small. Recently, the decay rate of the SM-like vacuum has been studied in 
detail for such a case [IHl 133 SSI SS] . On the other hand, if /i is as large as the stop masses, 
the down-type Higgs boson also has a large VEV at the CCB vacua due to the tri-linear 
scalar coupling among stops and the down-type Higgs, which is proportional to ytfj.. The 
vacuum stability condition is important in such a case because signihcant deviations of the 
Higgs partial widths from the SM prediction may occur. In order to study the SM-like 
vacuum stability, we consider the tree-level scalar potential in the held space involving ti, 
t^, hu and (which are canonically normalized scalar helds embedded in the left-handed 
stop, right-handed stop, the up-type Higgs and the down-type Higgs, respectively). 

The relevant part of the potential is given by 


V =^m^n hi + hi - m\^ ^r i% 

+ -^yMthu - fihd)iLiR + + Pihl + ilhl) 


+ y/M - Pgf + YsiCi 


hl-pL? + Y2Sy{hl-hl + lt 
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where 


(2.25) 


ml^ = m\ sin^ j3 — 


( 2 . 26 ) 


(2.27) 

(2.28) 


m22 = rn\ cos^ /3 + 

X 

"^12 = sin 2/3. 


When the CCB vacua become deeper than the SM-like vacuum, the latter is not stable. 
The vacuum decay rate per unit volume is calculated in the semi-classical approximation, 
and is expressed as 


r/1/= C'exp(-5i5), (2.29) 

where S'e is the Euclidean action of the so-called bounce solution [SHl For the decay 
of the SM-like vacuum, C is estimated to be ~ (100 GeV)^. In order for the lifetime of the 
SM-like vacuum to be longer than the present age of the universe, the Euclidean action is 
required to satisfy 


Se > 400. 


(2.30) 


In our numerical analysis, CosmoTransition 2.0al [52] is used to hnd the bounce solution 
in the four-dimensional held space parameterized by hd, hu, ti and fj?, and to calculate the 
Euclidean bounce action Se- The calculation is done at zero temperature, and therefore 
thermal effects are not taken into account. The model parameters such as the tri-linear and 
top Yukawa couplings are evaluated at the SUSY scale Msusy- 

In Fig. |2]we show the contours of constant Se for the PS and PL solutions on ruA vs. /i 
plane. Here, tuq = rrijy = M 3 = 5TeV, while sfermion masses, mi and m^, are taken 
to be max(mp,/i) (see Sec. |3|). In each plot, the region between the solid (dashed) green 
lines satishes Se > 400 (300). For the PS solution, these lines appear only for |/i| > lOTeV. 
On the other hand, the upper bound on |/i| is comparable to the stop masses for the PL 
solution. 

The tree-level potential Eq. fl2.25p is used for our numerical analysis. Radiative correc¬ 
tions may change the scalar potential particularly around the SM-like vacuum. In order to 
study their effects, we also estimated the decay rate of the SM-like vacuum by including 
the top-stop and bottom-sbottom one-loop correction to the Higgs potential]^ We found 
that the Euclidean action tends to increase by about 10 - 20% for Se ^ 400 from that of 
the tree-level potential. In order to see the sensitivity of Rg to /r, we show the contours of 
Se = 300 in the same plot. 

When both /i and tan/3 are large, we might have to consider other CCB vacua in the 
sbottom-Higgs direction, which is driven by the tri-linear coupling of yb{Abhd — iJ,hu)bLbji /\/2 
(where b^ and bji are left- and right-handed sbottoms, respectively) (cf. Ref. [H]). If both 
fi and Ab are large, the bottom Yukawa coupling can be enhanced not only by tan/3 but 
also by (1 -|- A;,)”^ (see Eq. fl2.3ip L However, when the squark masses are universal, such 


A better treatment would be to introduce full one-loop radiative corrections to the potential involving 
stop, sbottom and the Higgs boson so that the potential is stable against the renormalization scale at least 
at the one-loop level. 
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Figure 2: Contours of Se = 300 and 400 on vs. jj, plane for the PS and PL solutions of 
At. Here, tuq = = 5TeV, rrij^ = = max(mf/, |/r|), and tan/? = 20. 


a parameter region is already excluded by the other constraints discussed in this section. 
Therefore we do not consider the constraint coming from the CCB vacuum involving the 
sbottom sector in this paper. 


2.2.3 Bottom Yukawa coupling 


When tan/3 is very large, the bottom Yukawa coupling is sizable. We can impose an up¬ 
per bound on tan/3 by requiring perturbativity of the Yukawa coupling constants up to, 
for instance, the GUT scale. In the MSSM, the bottom Yukawa coupling constant i/b is 
proportional to (1 -|- and hence, yb is enhanced when Ab is negative (see Eq. fl2.4p i. 

Consequently, the bound on tan /3 is more severe when Ab < 0. 

In our numerical analysis, we estimate the bottom Yukawa coupling constant using the 
following relation 


#7 


//b(MsuSY) ^ 


\/2mj,(MsusY) 
V cos /3(1 -|- Ab) 


(2.31) 


Then, we follow the evolutions of coupling constants by solving the renormalization group 
equations at the one-loop level and impose a condition that the bottom Yukawa coupling is 
perturbative up to the grand unihed theory (GUT) scale Mqet- Numerically, we require 


|?/6 (^gut)| < 1, (2.32) 

discussing the perturbative bound, we neglect holomorphic corrections to the Yukawa coupling con¬ 
stant, which are orders-of-magnitude smaller than the tree-level value of the Yukawa coupling constant for 

tand ^ 1. 
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where we take Mqut = 2 x 10^® GeV. This constraint excludes a large tan /3 region especially 
when Af, < 0. Notice that A;, is a non-decoupling parameter, and hence, this constraint is 
important even in the limit of heavy superparticles. 


3 Higgs Partial Decay Widths 

In this section,we discuss the partial decay width of the lightest Higgs boson h. We dehne 
the ratio of the partial decay width of the lightest Higgs boson to that of the SM prediction: 

T(h^F) 

where F denotes a specihc hnal state. In the following, we show how much Rp can deviate 
from the SM prediction [Rp = 1) for various hnal states. As we have mentioned, FeynHiggs 
2.10.2 is used to calculate the partial decay widths of the Higgs boson, in which full one- 
loop contributions are taken into account for the fermionic hnal states. In the package, a 
resummation of the Ah corrections is also included in calculating the partial decay width for 
h^bb^. 

In our numerical calculation, we adopt (approximate) GUT relation among the SU{3)c, 
SU{2 )l, and U{1)y gaugino masses: 

Af3(MsusY) = 3M2 (Msusy) = 6 Mi(Msusy)- (3.2) 


All the phases in the MSSM parameters are assumed to be negligible, and we adopt the 
convention of M 3 > 0. For simplicity, we also assume that the sfermion masses are universal 
with respect to the generation indices. 

First, we show the soft mass dependence of Rp without taking the phenomenological 
constraints into account. In Fig. [3l Rp are shown for F = 77 , W^W~, gg, bb, and rr as 
functions oUmsoft, taking M 3 = rriA = —/^ = mg = mfj = mj^ = mi = = m^oit and 


tan (3 = 40|^ By choosing negative /r with \fi\ ~ rusoft, the correction to the bottom Yukawa 


coupling becomes signihcant. 

Remarkably, even though rusoft is around several TeV, the partial decay widths for F = bb, 
TT and gg deviate from the SM prediction by more than (9(1) % for the PL solution. However, 
SUSY contributions to the other widths are smaller. Note that Rih — 1 is as about twice 
as Rfr — 1 for the PL solution, which is due to the difference between Ah and A,-. We 
also note here that the partial decay widths have an appropriate decoupling behaviour, i.e., 
P(h —)■ F) ^ F) as the masses of superparticles and the heavy Higgs bosons 

become inhnitely large. 

Next, let us include the phenomenological constraints discussed in Sec. 12.21 In Fig. 01 
contours of Rih — 1 are shown on vs. tan/3 plane. Here, all the sfermion masses, M3, 
and |/i| are set to be 5 TeV. In the cases of the PS and PL solutions, the contours end at 


We have checked that ~ Rzz holds in the parameter region of our study. 
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msoft [GeV] msoft [GeV] 


Figure 3: Ryv (left) and Rjj (right) are shown as functions of rusoft; here all the fermion 
masses, M 3 , and |/x| are equal to msoft, and tan /3 = 40. (The sign of fi is taken to be negative, 
while Wino and Bino masses are given by using the approximate GUT relation fl3.2p .i In 
the left (right) plot, the black, red, and blue lines correspond to h —)■ 77 ( 66 ), h —)■ W^W~ 
(fr), and h —)■ gg, respectively. The solid (dashed) lines correspond to the PL (PS) solutions 
of At. Note that the constraints in Sec. 12.21 are not taken into account. 

around tan/? ~ 45. This is because, when tan/? is large, there is no solution for At to satisfy 
ruh = 125.7 GeV. The bottom-sbottom loop contribution interferes destructively with the 
top-stop contribution, and thus, larger (smaller) At is required for the PS (PL) solution. 
Then, and Af^'^ merge into a single solution at certain tan/?, which is the value where 
the top-stop contributions cannot raise the Higgs mass anymore for hxed msoft- Note that, 
Afe is negative when /r < 0 and At > 0. In such a case, yb is signihcantly larger than the 
tree-level value, and thus, the bottom-sbottom loop contributions are enhanced. 

In Fig. m the constraints from Br(i?s —)■ gAy~), the vacuum stability, and the perturba- 
tivity of the bottom Yukawa couplings are shown. Wide parameter region is excluded when 
At is large, i.e., in the PL and NL panels. The constraints from Br(i?s —)■ gAgr) and the 
vacuum stability become weaker for large m^, while that from the perturbativity is not. 
For the PL solution, the constraints from Br(i?s —)■ g'^g~) and the vacuum stability change 
drastically as we vary tan /? for tan /? > 30, where the bottom Yukawa coupling is much 
larger than the tree-level value. 

Even if the masses of superparticles are relatively large (i.e., 5 TeV), —1 can be as large 
as 0{1) %. Such a large deviation may be within the reach of future collider experiments. 
Expected accuracies at the future experiments have been discussed (see Table 1-16 of Ref. [3] 
and Ref. H); the numbers are summarized in Tabled] The accuracies of 5T[h —)■ 66 ) = 
0.64%, 5V{h —)■ fr) = 1.1%, and 5V{h —)■ gg) = 1.2% are claimed to be achievable at e+e“ 
colliders ultimately. Therefore, it is found that, even if the superparticles are kinematically 
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Figure 4: Contours of Rii, — 1 are shown for the PS, NS, PL, and NL solutions of At. Here, 
all the sfermion masses, M 3 , and |/i| are taken to be 5TeV, and the sign of /i is set to be 
negative. The left regions of the blue lines are excluded by Bt{Bs —)■ while those of 

the green solid (dashed) lines are constrained by the vacuum stability condition, Se > 400 
(360). The bottom Yukawa coupling becomes non-perturbative below the GUT scale in the 
region above the red line. 
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LHG[3] 

iLcm 

TLEP[3] 

^/s [GeV] 

1400 

1400 

250 

500 

1000 

240 

350 

fdtC[fh ^ 


300 

3000 

250 

500 

1000 

2500 

10000 

+2600 

77 


10 - 14 

4-10 

38 

17 

5.8 

3.8 

3.4 

3.0 

99 


12 - 16 

6-10 

12 

4.0 

1.6 

1.2 

2.2 

1.6 

66 


20 - 26 

8-14 

9.4 

1.9 

0.78 

0.64 

1.8 

0.84 

fr 


12 - 16 

4-10 

10 

3.8 

1.6 

1.3 

1.9 

1.1 


Table 1: Expected accuracies of the determinations of the partial decay widths of the Higgs 
boson in units of percents mm- The accuracies of the widths are assumed to be twice the 
accuracies of the determinations of couplings. 


unaccessible at the LHC, we may observe the MSSM signal by studying the partial decay 
widths of h in detail. 

For the PL solution, it is also found that the partial width of h —)■ 66 can deviate from the 
SM prediction by about 2% even for ttia = 6 TeV. On the other hand, we checked that Rfr 
is smaller by about 1 % than at the same parameter point, because the non-holomorphic 
correction |Af,| is bigger than |At-| at this model point. When tan/3 is smaller, the difference 
between and Rfr decreases, since A^ is approximately proportional to tan/3; in such a 
case both R^, and Rfr are well approximated by (sina/ cos/3)^. 

Let us see how much Ri^, and Rfr can change in the parameter space consistent with the 
phenomenological constraints. We have performed the scan in the following parameter space 
of the MSSM: 


rriQ = rrijj = M 3 = 2, 3, 4, and 5 TeV, 
rrif) = mi = = max(mp, |/i|), 

A _ 4 (NS) .(NL) .(PS) .(PL) 

1 1 1 1 


• 0.8 TeV < mA < 6 TeV, 

• —5 < fi/mjj < —0.5, or 0.5 < /i/mp < 5, 

• 5 < tan/3 < 50. 

Here, At is determined to satisfy mh = 125.7 GeV. The other tri-linear couplings (e.g., Ab 
and Ar) are assumed to be equal to At, we checked that our numerical results are insensitive 
to this assumption. We take = mp = M 3 , while sfermion masses other than rUg and mp 
are set to be equal to max(mp, |/i|). We have checked that Rib and Rfr are almost insensitive 
to these scalar masses unless they are very small. However, when |/r| is much larger than mp, 
mp, and mp, bottom and stau mixings become sizable which causes additional complexity 
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Figure 5: The ranges of R^f, (blue) and Rfr (red) with small \At\ solutions (i.e, the PS and 
NS solutions for At). We take ttIq = rrijj = M 3 = 2, 3, 4 and 5TeV. Other sfermion 
masses are set to be max(mQ, |/i|). The partial widths are constrained by Br(i?s —)■ /i+/i“), 
the vacuum stability condition {Se > 400), and the perturbativity of the bottom Yukawa 
coupling (|?/ft(GUT)| < 1). The black region corresponds to the tree-level prediction for 
tan/3 = 5-50. The black dashed line is the SM prediction. 
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Figure 6: Same as Fig. |5l but with all the solutions of At. 
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to our parameter scan, and/or lighter stau become tachyonic. The values of mi, and 
are chosen to avoid this problem. 

If \fi\ would become much larger than squark masses (e.g., |/i|/mp S> 5), the Higgs partial 
widths could deviate from the SM prediction sizably, since the Higgs mixing angle would 
be enhanced by radiative corrections. However, the vacuum stability condition excludes 
signihcant amount of the parameter region with large |/i|, as we have shown in Fig. |2J In 
our scan, we set |/r|/mp < 5; this upper bound is large enough to cover the whole region 
allowed by the vacuum stability. 

In Figs. |5] and El we show the minimal and maximal values of Rii, and Rfr as functions of 
mA and rrisoft, taking into account the phenomenological constraints discussed in the previous 
section. Fig. |5] corresponds to the small At solutions (NS and PS). The largest and smallest 
values of R^, and Rfr are achieved when \fi\ and tan/5 are large and are marginally allowed 
by the phenomenological constraints. The deviations mainly come from radiative corrections 
to the Higgs mixing angle and Ab- In particular, the maximal value of Rib is larger than 
that of Rfr for large mg = mjj. As discussed earlier, Ab becomes negative and sizable with 
fi<0, which results in a signihcant enhancement of Rib- 

Fig. El shows the minimal and maximal values of Rib and Rfr, taking all the solutions of 
NS, PS, NL, and PL into consideration. They change signihcantly compared to Fig. El The 
partial widths become extremum when At takes the NL or PL solution except the maximal 
value for small mA- For < 1 TeV, the NS or PS solutions give the maximal value, since 
the phenomenological constraints, especially the vacuum stability condition, are too severe 
for the NL and PL solutions (see Fig. 0]). As in the case of Fig. El the maximal value of 
Rib is much larger than that of Rfr because lA^I S> |At-|. The maximal value of Rib has 
a non-trivial bump-like structure. When mA is small, medium (i.e., around the peak of 
the bump), and large, the partial decay width is bounded by the vacuum stability, flavor 
constraint, and the perturbativity of yb, respectively. 

The measurement of the Higgs couplings may provide an evidence of BSM, because the 
values of Rib and Rfr may show signihcant deviation from the SM prediction. In future 
experiments, the partial decay widths of h ^bb and h ^ fr may be measured at the < 1% 
level (see Table E]). In the case of the PS and NS solutions. Rib — 1 can be as large as 2% 
(3%) if mA is 2.7TeV (2.2 TeV) for rusoft = 2 TeV, and 2.6 TeV (2.1 TeV) for rrisoft = 5 TeV. 
On the other hand, Rfr — 1 becomes larger than 2% (3%) when mA is smaller than 3.0 TeV 
(2.4 TeV) for rusoft = 2 TeV, and 2.1 TeV (1.8 TeV) for msoft = 5 TeV. Including the PL and 
NL solutions, deviations of 2 - 3% level are achieved with larger value of m^. For example, 
with TTisoft = 5 TeV, Rib — 1 = 3% can be achieved with mA = 6.0 TeV. 

In Fig. [3, we show the sfermion mass dependence of Rp with taking the phenomenological 
constraints into account. The parameter scan is performed in the same setup as Fig. El but 
mp is varied. Here, mA is hxed to be 5 TeV. In the left plot, the minimal and maximal 
values of Rib and Rfr are shown. They are achieved by the PL or NL solution. Since is 
hxed, the decoupling behaviour is not observed. Rather, the non-decoupling contribution to 
Af, enhances Rib as mp- increases, since larger \At\ is required to satisfy mh = 125.7GeV for 
the solutions. On the other hand, it is found that radiative corrections to the Higgs mixing 
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Figure 7: The ranges of Rjf (left) and Rgg (right) are shown as functions of rrify. The parame¬ 
ters are scanned in the same manner as Fig. [HI but tua is hxed to be 5 TeV while mp is varied. 
Here, all the solutions of At are included in the scan with imposing the phenomenological 
constraints discussed in Sec. Ea The black dashed line is the SM prediction. 


angle do not change so much even if mp increases. 

In the right panel of Fig. [TJ the minimal and maximal values of Rgg are displayed. The 
maximal value occurs for the PS or NS solution, while the minimal value is achieved by the 
PL or NL solution. We can see the decoupling behavior, i.e., Rgg approaches to unity as 
mp increases. The partial decay width of h ^ gg can deviate from the SM prediction by 
about 3% for mp = 2 TeV, while it decreases rapidly and becomes 1% for mp = 3.7 TeV. 
We have also checked that these values do not change so much for m^ = 8 TeV. (However, 
they change signihcantly if niA is smaller, since the phenomenological constraints exclude 
the parameter space severely.) According to Table [H the partial width of h ^ gg is expected 
to be measured at the 1.2% accuracy in future experiments. Thus, as far as superparticles 
are relatively light, we may observe a signal of the MSSM in the measurements of this partial 
width even if the heavier Higgses are out of the reach of the LHC. 

Finally, we show how large fraction of the parameter space can be covered by future e’''e“ 
colliders. For this purpose, we dehne the 6xf variable as 


^ 2 \T{h^ F)F)? 

~ [5T{h F)]2 


(3.3) 


where (5F(h —)■ F) is the expected accuracies of the determinations of the Higgs partial 
decay widths at ILC with ^/s = 1 TeV and JdtC = 2500 fb“^ (see Table E]). Based on 
this quantity, we dehne the parameter region which is accessible with ILC at Sxp > 4. We 
perform a parameter scan and study if each model point is accessible with ILC and satishes 
the phenomenological constraints. The MSSM parameters are scanned in the ranges of 
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Figure 8 : Upper-left: The number of model points accessible with ILC by at least one decay 
mode of h as a function of niA (green histogram), as well as that of model points allowed by 
the phenomenological constraints (dotted histogram). Upper-right: The number of model 
points allowed by the phenomenological constraints on niA vs. tan/? plane. Lower-left: The 
number of model points accessible with ILC by h — bb. Lower-right: The number of model 
points accessible with ILC by h —rr. 


800 GeV < niA < 5 TeV (with the step size of 100 GeV), 5 < tan/3 < 50 (with the step size 
of 1), and 0.5 < |/i|/mp < 2 (with the step size of 0.1). Thus, for each set of (m^,tan/9), 
384 model points are studied, taking account of positive and negative values of /i as well 
as all four solutions of At. In addition, we adopt the relations mg = M 3 = mp, and 
mp, = mp = mp = max(mp, |/i|). Goncentrating on the parameter space where the LHG 
will have a difficulty in hnding superparticles, mp is taken to be 3, 4, and 5 TeV. In Fig. [HI 
the ILG coverage of the parameter space of our scan is displayed. In the upper-left panel, we 
show the distribution of the number of model points accessible with ILG by any of the decay 
modes of h as a function of rriA. The green histogram is a distribution of the number of 
model points which satisfy 6xf > 4 and the phenomenological constraints, while the dotted 
one is that satisfying the phenomenological constraints without imposing Sxj? > 4. We hnd 
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that the number reduces drastically at tua ~ 2 TeV. In the upper-right panel, we show the 
number of model points which survive the phenomenological constraints on vs. tan/? 
plane. Here, 5x1^ > 4 is not imposed. Then, in the lower panels, we show the numbers 
of model points which can be accessed by the decay modes of h ^ bb (lower-left) and 
h —)■ fr (lower-right) with taking account of the phenomenological constraints. They reduce 
signihcantly for > 2 TeV irrespective of tan f3. The accessible points for tua ^ 2 TeV are 
mostly with PL or NL solution. 

4 Conclusions and Discussion 

In this paper, we have studied the partial decay widths of the lightest Higgs boson in the 
MSSM. Taking account of relevant phenomenological constraints, i.e., Higgs mass, flavor 
constraints, vacuum stability, and the perturbativity of coupling constants up to the GUT 
scale, we have calculated the expected deviations of the partial decay widths from the SM 
predictions. 

The partial decay widths are enhanced if the /i-parameter is relatively large. However, 
such a choice may conflict with some of the phenomenological constraints. In particular, the 
vacuum-stability condition imposes a stringent constraint on the parameter space. We have 
found that, with too large \fj,\, there show up CCB vacua where the down-type Higgs held as 
well as the up-type Higgs and stop helds acquire large VEVs; existence of such CCB vacua 
was not seriously considered in the previous studies. In addition, when /rtau/S is large, 
non-holomorphic correction to the bottom Yukawa interaction becomes so large that the 
bottom Yukawa coupling constant becomes non-perturbative below the GUT scale. Large 
value of /itau/S may also cause too large havor-violating decay of H-mesons. By taking 
these constraints into account, the maximal and minimal possible values of the Higgs partial 
widths are restricted. 

We found that the deviations of the partial decay widths from the SM predictions can 
be of 0{1) % for some of the decay modes. In particular, those of T{h —)■ bb) and T{h —)■ fr) 
may show signihcant deviations even if the superparticles are out of the reach of 14TeV 
LHC. In addition, the deviation of r{h —)■ gg) may also be sizable if the superparticles are 
relatively light. We emphasize that, although our scan is limited to some part of the MSSM 
parameter space, we have found the regions where the deviations from the SM predictions 
are within the reach of proposed e’''e“ colliders even if superparticles would not be observed 
at the LHC. 
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